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Abstract. We study the linear, but fully non-adiabatic tidal re- 
sponse of a uniformly rotating, somewhat evolved (X c = 0.4), 
10 Mq main sequence star to the dominant I = 2 components 
of its binary companion's tidal potential. This is done numer- 
ically with a 2D implicit finite difference scheme. We assume 
the spin vector of the 10 Mq star to be aligned perpendicular 
to the orbital plane and calculate the frequency a and width 
of the resonances with the prograde and retrograde gravity (g) 
modes as well as the resonances with quasi-toroidal rotational 
(r) modes for varying rotation rates Sl s of the main sequence 
star. For all applied forcing frequencies we determine the rate 
of tidal energy and angular momentum exchange with the com- 
panion. In a rotating star tidal energy is transferred from I = 2 
g-modes to g-modes of higher spherical degree (I = 4, 6, 8, . . .) 
by the Coriolis force. These latter modes have shorter wave- 
length and are damped more heavily, so that the I — 2 resonant 
tidal interaction tends to be reduced for large rotation rates f2 s . 
On the other hand, the density of potential resonances (a broad 
I spectrum) increases. We find several inertially excited unsta- 
ble I > 4 g-modes, but not more than one (retrograde) unstable 
I = 2 g-mode and that only for rapid rotation. Our numerical 
results can be applied to study the tidal evolution of eccentric 
binaries containing early type B-star components. 



tide. To this end a 2D implicit code was developed for which 
the effects of the Coriolis force on the non-radial oscillations 
are taken fully into account. These earlier studies of rotational 
effects were based on a chemically homogeneous 20 M© star 
and were focussed on the low-frequency inertial regime, i.e. to 
forcing frequencies a < 251 s , where a is the forcing frequency 
in the frame corotating with the star at a rate tt s . The present 
study extends these calculations to a somewhat evolved (with 
core hydrogen abundance X c = 0.4) 10 M star whereby the 
forcing frequency a runs from low (high radial order g-modes) 
to high frequencies (up to gi). The motivation for this work is 
that we intend to apply the calculated energy and angular mo- 
mentum exchange rates with the companion to study the tidal 
evolution of eccentric early type binary systems. Recently sev- 
eral new studies of tidal evolution in eccentric early type bi- 
naries (with a com pact comp anion) have appeared (e.g. Kumar 
& Goodman 1996; Lai 1997). However, these studies have not 



considered the important effect of rotational (r) modes on the 
tidal exchange of angular momentum in eccentric binary sys- 
tems. We anticipate interesting tidal effects in eccentric bina- 
ries by the counteracting effects of resonant prograde g-modes 
and retrograde r-modes when the early type component is rotat- 
ing near its 'pseudo' synchronous rate at periastron. This will 
be studied in a following paper. 
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1. Introduction 



Stars: binaries: close- 



Zahn ( 1977 ) initiated the study of radiative damping of the dy- 
namical tide as a viable mechanism for effective tidal interac- 
tionin early type close binary systems. Savonije & Papaloizou 
( 1984 ) were the first to perform fully non-adiabatic calculations 
of the dynamical tide and to study the interplay between stellar- 
and tidal evolution which appears crucial in understanding the 
effects of tid es in early type stars. More recentl y, Pap aloizou 
& Savonije ( 1997 ) and Savonije & Papaloizou ( 1997 ) (from 
now on SP97) studied the effects of rotation on the dynamical 
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2. Basic equations 

We consider a uniformly rotating, early type star with mass M s 
and radius R s in a close binary with circular orbit with angular 
velocity u> and orbital separation a. We assume the stellar an- 
gular velocity of rotation fl s to be much smaller than its break- 
up speed, i.e. (f2 s /O c ) 2 < 1, with SI 2 = GM S /R^, so that 
effects of centrifugal distortion (oc SI 2 ) may be neglected in 
first approximation. We wish to study the response of this uni- 
formly rotating star to a perturbing time-dependent tidal force. 
The Coriolis acceleration is proportional to Sl s and we take its 
effect on the tidally induced motions in the star fully into ac- 
count. We use spherical coordinates (r, ip), with origin at the 
stellar centre, whereby $ = corresponds to its rotation axis 
which we assume to be parallel to the orbital angular momen- 
tum vector. We take the coordinates to be non-rotating. 
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As is well known, in a non-rotating star the solutions 
of the linearized non-radial stellar oscillation equations can 
be expressed in terms of spherical harmonics, i.e. the spa- 
tial part of each mode can be fully separated into r-, 
and (^-factors (e.g. ^edoux & Walraven 1958| ) U (r, d, tp) = 



u(r)P z m (cos tf) e" 



where P, m represents the associated 



Legendre polynomials for I and m. 

The introduction of the Coriolis force, however, destroys 
the full separability of the oscillation equations, it only being 
retained for the (/^-dependence. It turns out (e.g. Berthomieu 
et al. 1978) that two independent sets of approximately 
spheroidal oscillation modes exist: modes in which the density 
perturbation is even with respect to reflection in the equatorial 
plane, which have I — \m\ even valued, and modes with odd 
symmetry for the density, having I — \m\ odd valued. In ad- 
dition, for each /, there is a set of quasi-toroidal r-modes (e.g. 
Papaloizou & Pringle 1978 ) which couple with the spheroidal 
modes of I ± 1. 

Let us denote perturbed Eulerian quantities like pressure 
P', density p', temperature T 1 and energy flux F 1 with a prime. 
The linearized hydrodynamic equations governing the non- 
adiabatic response of the uniformly rotating star to the perturb- 
ing potential $t may then be written 
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where e, are the unit vectors of our spherical coordinate sys- 
tem, k is the unit vector along the rotation axis, v' denotes the 
velocity perturbation, k the opacity of stellar material and S 
its specific entropy. These perturbation equations represent, re- 
spectively, conservation of momentum, conservation of mass 
and conservation of energy, while the last equation describes 
the radiative diffusion of the perturbed energy flux. For sim- 
plicity we adopt the Cowling ( 1941 ) approximation, i.e. we ne- 
glect perturbations to the gravitational potential caused by the 
stellar distortion. We also neglect perturbations of the nuclear 
energy sources and of convection. 

For a circular orbit (with orbital angular speed u>) the com- 
panion's perturbing potential can be expanded as the real part 
of (e.g. |Morse & Feshbach 1952[): 
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where M p is the companion's mass, a the orbital separation, 
P ; m (cos d) the associated Legendre polynomial and e m = 1 
for m = and 2 for m > 0. We will consider only the domi- 
nant I = 2 components of the tidal forcing. Adopting the same 
azimuthal m symmetry and time dependence for the perturbed 
quantities as the forcing potential, the perturbed velocity can 
be expressed as v' = ia£, where a = m((J — fl s ) is the forcing 
frequency felt by a mass element in the uniformly rotating star 
and £ is the displacement vector. 

The perturbations can be written as e.g. £ r (r, p, t) — 
£ r (f) e lm( - wt_¥ '- ) where ^ r is the radial component of the dis- 
placement vector, while £ r (?", i?) is assumed complex to de- 
scribe the azimuthal phase shift with respect to the forcing po- 
tential (||) induced by any occurring dissipation, e.g. turbulent 
viscosity or radiative damping, see energy equation below. 

The current equations contain extra terms compared to 
those in SP97 because of the occurring mean-molecular weight 
(/^-gradients near the edge of the convective core. We assume 
diffusive mixing to be negligible on the (oscillation) timescales 
under consideration, so that the Lagrangian variation of the 
mean molecular weight Sfi a — or ^ = — dl ^ a £ r . We can 
thus use 
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to eliminate the pressure perturbation, where \p = g 1 ;" P . 

Xt — g\nT anc ' = g/"^ follow from the equation of 
state. 

Writing for simplicity from now on £ r for £r(fj#)> etc., 
while dividing out the factor e im (" t_v ), Eqs. yield the 
seven scalar Eqs. |7|-|l3| given below. 

First of all we write out the perturbed equation of continuity 
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We can use (J7|) to eliminate the term introduced by the 
radial derivative of the pressure perturbation (through (||)) from 
the radial equation of motion, so that the latter equation can 
be expressed (after adding viscous terms to introduce turbulent 
damping in convective regions, see SP97) as 
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where // = cos i3 and £ is the coefficient of turbulent viscosity 
defined below. The t?-equation of motion becomes 
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By applying the thermodynamic relation 
P 1 



SS = S' + £ ■ VS = 
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where the symbol 5 denotes a Lagrangian perturbation and Ta 
the adiabatic exponents of Chandrasekhar, the perturbed energy 
equation can be expressed as 
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where rj — (T3 — 1) -gp is a local characteristic radiative diffu- 
sion length in the star, with F the unperturbed (radial) energy 
flux. We have eliminated F' v with help of the (^-component of 
the radiative flux equation. In the stellar interior 77 ~ which 



corresponds to almost adiabatic response. However, even for 
'high' frequencies a/il c rj 1 the diffusion length 77 becomes 
comparable to the scale height when the stellar surface is ap- 
proached. The associated radiative energy losses give rise to 
damping of the tidally excited oscillations whereby the result- 
ing phase lag with the companion generates a torque. 
The perturbed radial radiative energy flux is given by 
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where n„ = 



the hydrogen abundance. 

Finally the ^-component of the perturbed flux follows as 
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2.1. Boundary conditions 



The differential equations are supplemented by the following 
boundary conditions: at the stellar centre we require £ r and F' r 
to vanish, while at the stellar surface, we require the Lagrangian 
pressure perturbations to vanish 



P[_ 
~~P 



dlnP 







(14) 



and the temperature and flux perturbations to fulfill Stefan- 
Boltzmann's law 



(15) 



Furthermore, and F'g must vanish on the rotation axis while 
in view of the symmetry of the tidal force we adopt mirror sym- 
metry about the equatorial plane, i.e., for 1) = ir/2 we also 
require £0 and F'g to vanish. 

2.2. The unperturbed stellar model 



A recent version (Pols et al. 1995) of the stellar evolution code 



developed by Eggleton (1972) was used to construct the un- 
perturbed stellar input model. The model represents a some- 
what evolved main-sequence star of 10 M© with core hy- 
drogen abundance of X = 0.4 and Z = 0.02. The mass 
in the convective core is approximately 2.1 Mq. The model 
comprises 1200 (radial) zones and was constructed with the 
OPAL opacities (Iglesias & Rogers 1996). The stellar radius 
equals R s 
ature T e g 



— 3.825 x 10 TT cm, while the effective temper- 
= 2.314 x 10 4 K and the stellar moment of in- 



ertia I s = 1.56 x 10 56 gem 2 . The break-up angular speed 
equals £l c = 1.54 x 10~ 4 s _1 . The Brunt-Vaisala' frequency 

A = JiF (pHf _ TTF^) 8 iven in units of n ^ is P lotted in 
Fig. [y. The Brunt-Vaisala frequency attains large values in the 
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Fig. 1. Characteristics of the stellar model: continuous curve 
represents the Brunt-Vaisala frequency ^bv = sign^4-^/|^4| in 
units of the stellar break-up speedO c as a function of radial 
mesh number. The dashed curve shows the forcing frequency a 
for which the corresponding oscillation period equals the local 
thermal timescale. 



region where the convective core has retreated during the evo- 
lution and a composition gradient is formed (the 1 p a -gradient 
zone'). It can also be seen that there are two shallow convec- 
tive shells near the stellar surface. The dashed curve shows the 
'thermal frequency' v t h — 27r/r t h, where r t h = ^^Em ~ 
is a local thermal timescale, n is the opacity, c the velocity of 
light, and (3 the ratio of gas to total pressure. Non-adiabatic ef- 
fects become important at locations where f t h is comparable or 
larger than the forcing frequency a. It can be seen in Fig. [j] that 
this is the case for the layers in and above the inner convective 
shell. At the inner edge of this shell i/ t j, ~ Vl c , at the outer edge 
v th ~ 50 c . 

2.3. Turbulent viscosity in convective regions 

In convective regions the equations of motion are supplemented 
by extra terms to account for the occurring turbulent dissi- 
pation, as described in PS97. For the coefficient of turbulent 
viscosity ( we adopt a simple local mixing length approxi- 
mation: ( = a Hp v c , where H p is the pressure scale height, 

a = 2 is the mixing length parameter, and v c = ( ^ 
a characteristic convective velocity. For high oscillation fre- 
quencies the viscosity is reduced ( poldreich & Keeley 1977 ) 

by a factor min (l, (<tt c ) where r c = aH p /v c is the 
convective timescale. We limit £ to be everywhere less than 
5 x 10 12 cm 2 s _1 . 

In the /j a -gradient zone adjacent to the convective core, 
where the Brunt-Vaisala fequency attains large values, the tidal 
response has short wavelength and cannot be resolved on the 
grid when the forcing frequency \a\ drops to small values. We 
retain some viscosity in this region by letting the viscosity de- 
cay outwards from the boundary of the convective core (at r c ) 



as £ oc exp — [(r — r c )/(0.1 H p )] 2 . Note in this respect that, 
although the composition gradient suppresses radial overshoot- 
ing, horizontal turbulence may be well developed in this bound- 
ary layer between the core and envelope. Also the oscillation 
amplitudes are relatively large in this region. 

3. The torque integral: transfer of energy and angular 
momentum 

For a circular orbit the tidal potential (^) has nom = compo- 
nent. However, since we wish to apply our results to eccentric 
binaries for which the tide has an axisymmetric time dependent 
component we replace the factor e lm ( wt -¥ > ) by e 1 ( crt - m v) so 
that we can study m = forcing. Once we have solved Eqs. 0- 
[l3| for a given stellar rotation rate fl s and a given (l,m,a), i.e. 
for a term 



-fi m r l Pr (cos #) cos(at- mip) 



(16) 



in the forcing potential (||), where /j m oc M p /a l+1 , the rate 
of angular momentum exchange with the companion's orbital 
motion can be calculated as an integral of the tidal force per 
unit volume Fi m = — pV$/ m over the volume of the star 
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= Re 



m p' (r, $) r 2 sin i9 di3 dip dr 



dp 



= rrnrfi m / / Im(p') (cos -d)r ' +2 sintfdtfdr 



= mT, n 



(17) 



where Rc and Im stand for the real and imaginary part. Since 
the torque and the stellar and orbital rotation vectors are all 
aligned, only the magnitude of the transferred spin and orbital 
angular momentum needs to be considered. The rate of tidal 
energy exchange with the (10 Mq) star can be expressed as 



ldV 



Re / p V$ /m v'dV-Re p'V<5> lm ■ v mt dV. (18) 



Substituting v' = i<r£ and v Iot = fl s r sin$ e v we find, after 
taking the real part, that the second integral on the right hand 
side equals mfl s Ti, n . The first integral on the right hand side 
can be rewritten by applying the equation of continuity (Q). Par- 
tial integration then shows that this term is also proportional to 
% m and equals a% m . Concluding, we find for the tidally in- 
duced rate of respectively the energy and angular momentum 
change in the star: 



E s — oTi m 



(19) 
(20) 



where the torque integral is given by 

T lm = nf im J ^ Im (p'(r, p)) P^fjj) r l+2 dp dr (21) 
with /i = cos i?. 
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Fig. 2. Torque integral % m versus forcing frequency a for forcing with I = 2 and m — 2 on a 10 M Q star rotating at twenty 
percent of breakup speed. Prograde and retrograde gj_ fc -modes with k = 4 and 10 are labeled, as well as the location of the 
r-modes. Crosses denote calculated points, while the drawn continuous curve represents a fit, see text. 



4. Numerical procedure 

The set of partial differential Eqs. |7|-p"3| with boundary condi- 
tions ( fl4| ) and ( |l5| ) is approximated by a set of finite difference 
equations (FDE's) on a 2D grid in r and i.e. on a meridional 
plane through the star. The FDE's are similar (except for the 
extra terms associated with the /i a -gradient) to those described 
in SP97 and are solved by the same implicit scheme as given 
in that paper. We use a grid of 1200 zones in the radial and 
46 zones in the d direction, analogous to the one used in SP97. 
Since the radial grid is staggered, with part of the unknown per- 
turbations defined on the odd and the rest on the even zones, the 
effective number of radial zones is 600. 

For given values of (I, m, a) in the applied forcing potential 
(16) and for a given uniform rotation rate f2 s of the star we 
can thus obtain the stellar response, i.e. the complex valued 

perturbations ($), ($) and 

on the 2D grid in the meridional plane of the perturbed star. By 
multiplying these values with the common factor e 1 ( <rf-mi P) we 
obtain the non-adiabatic stellar response to the prescribed tidal 
forcing throughout the star. 

Because the 2D implicit numerical scheme involves many 
(complex) matrix invertions the solution requires large com- 
puter facilities and it is essential to use an economic method 
for tracing the many resonances with the free stellar oscilla- 
tion modes required to study the dynamical tide in a rotating 
star. The resonances are searched in frequency (<r) space by by 
tracing the maximum of the function ^>(a) = J^i j ' £ij 
where £ is the displacement vector, a * denotes complex 
conjugation and i and j run over the radial and $ grid, re- 
spectively. The maxima of "J/ are searched by using a robust 
method based on cubic interpolation (NAG library routine) 
which requires the calculation of the first derivative = 

Eij (&J • %- + €ij • ■ Fortunately, the derivatives 

H and ^t|- can be obtained cheaply once we have obtained 



the solution of the perturbations etc. For by differentiat- 
ing Eqs. |7j-|l3| with respect to a we arrive at the same equa- 
tions, except that the unknowns are now the derivatives of 
the perturbations (i.e. J| instead of £, etc.), while the right 
hand sides can be expressed in terms of the solved perturba- 
tions. The right hand side for equation (Q) for example becomes 
~2pa!;r — 2ipil s sin i? By storing the relevant inverted ma- 
trices during the solution procedure for the perturbations we 
then can simply combine the stored matrices with the adapted 
right hand sides to obtain a solution for SI at almost no extra 
costs. 

5. Results 

From now on we will express all frequencies in units of f2 c = 
1.54 x 10~ 4 s _1 . By varying the forcing frequency a = a — 
mO s ( m me corotating frame) for the I = 2 terms in the forcing 
potential (|6|) and by calculating the stellar response, we can 
search for resonances with the free stellar oscillation modes, 
applying the procedure described in the previous section. We 
adopted various values for the rotation rate of the 10 M star: 
fl s = 0, 0.1,0.2, 0.3 and 0.4. 

We traced the resonances with g^,-modes (up to k = 20) 
for both prograde and retrograde g 2 -modes. The lower index 
of g denotes k, the number of radial nodes in the displace- 
ment eigenvector, whereby the minus sign denotes a retrograde 
mode. Note that we characterize retrograde modes (i.e. modes 
that propagate in the direction counter to the stellar rotation) 
by negative values of the oscillation frequency a, i.e. we take 
m > for both prograde and retrograde modes. For rotating 
stars there is in addition a compact spectrum of r-modes, for 
slightly negative frequencies. The I = 2 component of the tidal 
potential excites r-modes with a predominant I = 3 component. 
Note that because we consider binary systems in which both 
stars are aligned perpendicular to the orbital plane the (/ = 2) 
forcing has only m = and |m| = 2 components, i.e. f 2 \ = 0. 
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Fig. 3. Torque integral 7j m versus (retrograde) forcing frequency a for forcing with ^ = 2 and to = 2 on a 10 Mq star rotating 
at fourty percent of breakup speed. Retrograde g 2 _ fe -modes with k = 4 and 10 (unstable) are labeled, as well as the r-modes with 
k = 2 and k = 5. Crosses and circles (unstable modes) denote calculated points. The peaks not fitted by the continuous curve 
correspond to resonances with higher spherical degrees I = 4, 6, . . . 



We can thus limit our calculations to either |to| = 2 or m = 
resonances. 

Note further that for all forcing frequencies a we keep 
the factor fi m in the forcing potential ( |l6| ) constant by adopt- 
ing a fixed value a = 4R S and M p = 1M©. Because we 
solve the linearized problem all listed values for the tidal 
torque integral ( plj ) can be simply scaled to the value for 
a required binary configuration by multiplying with a factor 
(M p /M Q ) 2 {m s /a) 2( - l+1 \ where a and M p are the actual or- 
bital separation and companion mass. Unless the companion is 
a compact star one should of course also take into account the 
(dynamical) tide in that star to determine the tidal evolution of 
the binary system. 



calculations are listed in Tables Al to A4 in the appendix. Note 



that for a non-rotating star the I = 2 spectrum is degenerate in 
to, so that for fl s = the resonance frequencies of I = to = 2 



g-modes are given by Table Al 



It can be seen that, as expected, for to = the g-mode res- 
onance frequencies do increase with il s , as is true for |oo| for 
the retrograde to = 2 g-modes and the r-modes. However, the 
prograde to = 2 g-mode resonances shift to lower oscillation 
frequencies as $7 S increases. Note in this respect that (approx- 
imate) conservation of radial vorticity, in combination with an 
increasing radial component of ambient rotation £! s cos d to- 
wards small co latitudes, tends to give retrograde wave propa- 
gation (see e.g. Unno et al. 1989). 



5.7. Resonance fitting 

Using the technique described in Sect. |J we determine the 
resonances with g- and r-modes and evaluate the torque in- 
tegral (|T]) from the calculated stellar response in the merid- 
ional plane of the 10 Mq star for all applied forcing frequen- 
cies. Close to each resonance we fit the obtained values for the 
torque integral by the resonance curve of a damped harmonic 
oscillator 



lm,0 



^ SoAs J 



(22) 



with eigenfrequency <7o, resonance full width Act and peak 
value Tim.Q. As a consequence of the Coriolis force, the stel- 
lar response to the I = 2 forcing contains a range of /-values 
(when we expand in spherical harmonics). However, we fit only 
the resonances with predominantly / = 2 (g-modes) or I = 3 
(r-modes). In between I = 2 resonances we assume % m (o) can 
be approximated by adding the contributions of the two adja- 
cent I = 2 resonances. The detailed results of our numerical 



5.2. Prograde and retrograde g-mode resonances 

In Figs. U and || we show the tidal torque integral % m versus 
forcing frequency a for a 10 M Q star rotating at a rate of 0.2 
and 0.4, respectively. Crosses and open circles correspond to 
numerical results, while the continuous curve corresponds to 
the above mentioned fit to individual resonances. When the 
10 M star rotates in the same sense as the orbital motion of its 
companion, the latter cannot excite modes with frequency (in 
the frame corotating with the star) less than a = -mf! s . Never- 
theless, we have calculated the whole range \a\ < 2. In Fig. ^|, 
the circles and the dashed fit curve represent unstable modes 
for which the torque integral has a reversed sign compared to a 



(see Sect. 5.4) 



The prograde g-modes shift to substantially larger oscil- 
lation frequencies compared to the unevolved ZAMS stellar 
model of the same mass. The retrograde g-modes shift to more 
negative frequencies, so that for ft s = 0.2 no strong retrograde 
g-modes can be tidally excited, unless the rather unlikely case 
of retrograde stellar rotation applies. However, when the stellar 
rotation rate is increased more and more retrograde g-modes 
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are shifted into the 'tidal window': up to g^ 13 for fi s = 0.3, 
and up to g 2 _ w for Q s = 0.4. 

But note that when the binary is eccentric and the early type 
star's spin is approximately synchronized at periastron, strong 
(retrograde) r-modes can be tidally excited even for relatively 
low stellar rotation rates. 

For small frequencies \a\ viscous damping in the fi a - 
gradient zone becomes significant as the local wavelength of 
the response gets short and the adopted turbulent viscosity gets 
large due to the com parable timescale of 'convection' and os- 
cillations (Sect. 2.3). Therefore the values for % m in the re- 



gion between g?_ 2 o an d me r-modes (and similarly for the cor- 
responding positive frequency range) depend on the uncertain 
assumptions about the viscous dissipation. Switching off all 
viscous dissipation in the ^-gradient zone yields for Q s = 0.4 
in the frequency range —0.4 < a < —0.2 torque values about 
an order of magnitude smaller than shown in Fig. H 



5.3. Effects of rapid rotation 

Although centrifugal effects are no longer negligible for the 
high rotation rates (O s = 0.4) considered here the results are 
still of interest, even though centrifugal distortion is neglected, 
because they show the effect of strong coupling by the Coriolis 
force. Note in this respect that in the higher density interior 
regions the local break up speed is substantially larger than fi c . 
A full calculation with centrifugal distortion included would 
require a much larger computing effort. 

The Coriolis force gives rise to strong coupling with modes 
of higher spherical degrees, so that for I = 2 forcing resonances 
with I = 4, 6, 8, . . . also become prominent, especially when 
located near a I — 2 resonance. Because the tidal forcing has 
I = 2 symmetry the latter resonances generally dominate, but 
all neighbouring resonances with / > 2 are also excited, so 
that the tidal response contains significant power over a broad 
range of (even) /-values. Often this makes mode identification 
(by means of a decomposition in Fourier-Legendre series of the 
eigenfunctions, see SP97) difficult since the power in different 
I -components are often comparable. One should consider the 
excited oscillation as a complex of coupled modes. Since, for a 
given frequency the g-modes with larger I have shorter wave- 
length, damping is enhanced, so that for large rotation rates the 
oscillation amplitude and the effective torque is reduced. This 
can be observed in Tables A1-A3, which show that the reso- 
nance area irAaTi, n of the I = 2 g-modes initially increases 
with il s but generally decreases when the rotation rate and in- 
ertial damping becomes large. 

The indirectly excited higher Z-resonances correspond to 
the peaks that have not been fitted in Fig. |3|. This figure shows 
part of the retrograde oscillation spectrum for il s = 0.4. Most 
of the stable unfitted peaks between g 2 L 3 and g 2 _ 5 correspond 
to what are (predominant) I = 4 resonances, except for the two 
nearest resonances around g 2 _ 3 , which are identified as g 8 _ 12 

and g 2 _ w are all 



and g_ 18 . The unfitted peaks between g_ ( 
(predominant) I — 4 resonances, and can be identified as gl 12 
to gl 18 consecutively. The g 2 _ & and g 4 L 14 resonances can hardly 



Table 1. Unstable m = 2 modes for f2 s = 0.4. The listed 
oscillation periods P m are relative to the inertial frame, a minus 
sign denoting retrograde propagation. 



mode Re(er) Im(cr) Pi n (d) 



-1.3826 -1.1 x 10 
-1.3451 -5.3 x 10 
-1.1524 -1.2 x 10 
-0.7602 -1.2 x 10 



-0.81 
-0.87 
-1.34 
11.86 



7T- 

4l:i 
4s 

it 8 
4s 



1.0580 -4.6 x 10~° 
1.1261 -4.3 x 10" 5 
1.4610 -1.0 x 10~ 4 
2.1887 -2.4 x 10" 5 



0.25 
0.25 
0.21 
0.16 



be distinguished from each other, they lie so close together that 
they form a strongly coupled complex, whereby the gl 14 reso- 
nance peaks one order of magnitude above the I = 2 resonance 
peak. Another symbiosis occurs for g 2 _ 10 and gl 18 , whereby, 
however, the former mode appears unstable and the latter sta- 
ble. In this frequency range the resonances with / > 4 are heav- 
ily damped as they approach their asymptotic low frequency 
regime. 

For still lower frequencies in the range around a = —0.5 
the calculated response in between resonances is significantly 
stronger than the wings of the fitted / = 2 resonance curves 
due to the smeared out contributions of higher I components. At 
these frequencies the spectrum of I — 4 modes starts to become 
dense whereby radiative and viscous damping gets heavy due 
to the short wavelength of the response. 

5.4. Unstable modes 



It is now known (e.g . Dziembowski & Pamyatnykh 1993 
|Gautschy & Saio 1993[ ) that most main sequence OB stars show 
unstable non-radial modes driven by the K-mechanism associ- 
ated with the metal opacity bump around 1.5 x 10 5 K in the 
OPAL opacity tables ([glesias & Rogers 1996). Two necessary 
conditions for a mode to be unstable by the k mechanism are 
that the Lagrangian pressure perturbation SP/P must reach a 
maximum value (coming from the stellar centre) in the driv- 
ing zone and that 2n/a is of order the thermal timescale in 
that region (e.g. Dziembowski et al. 1993). Obviously, when 
l^l <C vth = 27r/r t h in the ionization region that region is 
strongly non-adiabatic and no driving can occur. The 'thermal 
frequency' v t h (Sect. 2.2) in the lower convective shell asso- 



ciated with the opacity bump varies between approximately 1 
and 5. 

When the stellar rotation rate is increased some inertially 
excited (I > 2) modes in our 10 M model appear which have 
a torque integral % m with a sign opposite to that of er, which in- 
dicates instability. This can be checked by forcing the star with 



a complex frequency ( Papaloizou et al. 1997 ) and searching for 



maximum amplitude in complex frequency space. It appears 
that the fundamental p§ mode with |er| ~ 3.3 is unstable for all 
rotation rates Q s = to 0.4. There are more unstable p-modes 
which we have not studied. Table [j] lists the unstable g-modes 
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I = 2 



n, = 0.2 a. 



X r = 0.4 



41 
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-0.065 
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-0.055 



-0.05 



Fig. 4. Spectrum of r-modes: torque integral Ti m versus forcing 
frequency a for forcing with I = 2 and m = 2 on a 10 Mq star 
rotating at twenty percent of breakup speed. Crosses denote 
calculated points, the drawn continuous curve represents a fit. 



(prograde and retrograde) that are found for the highest rotation 
rate fl s = 0.4. 

Only for the highest rotation rate (f2 s = 0.4) an unsta- 
ble I = 2 g-mode appears: the retrograde gL 10 mode, with 
a — —0.76, which is located close to the stable gl 18 res- 
onance. The g 2 _ 10 mode has rather low frequency for the re- 
mechanism and is a bit peculiar in that it seems mixed with 
a strong short wavelength mode trapped in the \i a gradient 
zone. The oscillation period (in inertial frame) of the unstable 
I = 2 mode happens to be quite long: about 12 days, while the 
adopted stellar rotation period is 1.2 days. We also find four 
unstable prograde g-modes, see Table [jj. Our list of unstable 
modes is probably not exhaustive (except for I = 2), because 
we preferentially find those unstable modes which happen to be 
located near I — 2 resonances. Unstable modes, when excited 
by the tidal force, could give rise to tidal evolution counter to 
the normal direction of evolution. 



5.5. r-modes 

For small negative forcing frequencies around (e.g. Papaloizou 



& Pringle 1978): a 



2mO B 
1(1+1) 



strong resonances occur with 



quasi-toroidal oscillation modes analogous to Rossby modes 
in the earth's atmosfere (e.g. |Pedlosky 1979 ). For these modes 
the fluid elements oscillate almost exclusively in the horizontal 
direction whereby the restoring force is provided by the Corio- 
lis force (conservation of radial component of vorticity). Fig. Q 
shows the variation of the torque integral "7J m as the forcing 
frequency runs through the r-mode region for a star rotating 
with il s = 0.2. The adopted turbulent viscosity law in the fi a 
gradient zone adjacent to the convective core causes the r-mode 
resonances to become rather broadened compared to the results 
for ZAMS models (SP97). When we switch off the viscosity in 
the /i a -gradient zone the peak values for the torque become one 



or even two orders of magnitude higher, with a corresponding 
decrease of the resonance width. However, this will not make 
a large difference for the overall tidal evolution of the binary 
because the total area ~ 7rAa7j mi o of each resonance remains 
approximately constant. The resonant interaction with r-modes 
gets stronger with increasing rotation rate fi s . This can be seen 



in Table A.4 where the resonance area increases monotonically 
with fL. 



5.6. Excitation of inertial modes in the convective core 

Similar to the results obtained in SP97 for a 20 M Q ZAMS 
star, we find inertial wave oscillations in the convective core 
and in the convective shell region when the forcing frequency 
falls in the inertial range \a\ < 2il s . When we artificially cut 
the radiative envelope from our stellar model we find a dense 
spectrum of modes in the core. For some frequencies the reso- 
nances in the core appear very strong. The general amplitudes 
are an order of magnitude larger (or more) compared to the 
full stellar model where the inertial waves can leak out of the 
core. This would suggest that the inertial modes could possibly 
excite gravity waves in the envelope (SP97). However, when 
we consider a model in which the convective core is cut away, 
we find in this frequency region values for the torque integral 
comparable to the ones calculated for the full stellar model. 
Therefore it seems the inertial core modes do not significantly 
contribute to the stellar torque in this somewhat evolved stellar 
model. However, the oscillations in the /^-gradient zone out- 
side the convective core are poorly resolved for low frequen- 
cies, so that further studies with significantly more meshpoints 
in this region seem required to draw firm conclusions. 

6. Conclusions 

We have studied the linearized non-adiabatic tidal response of 
a somewhat evolved 10 M Q main sequence star to tidal forcing 
with spherical harmonics degree I = 2 and calculated the tidal 
exchange of energy and angular momentum with an orbiting 
companion as a function of stellar rotation rate taking the Cori- 
olis force fully into account. We found, as expected, that with 
increasing rotation rate the inertial coupling with higher spher- 
ical degree oscillation modes gets stronger, so that when the 
rotation rate becomes comparable to or larger than the tidal os- 
cillation frequency a significant amount of tidal energy is trans- 
ferred to these higher degree oscillations. Because of enhanced 
damping of the higher I oscillations fast rotation reduces the 
resonant tidal interaction with I = 2 g-modes. However, this 
is compensated by the appearance of an additional spectrum of 
(I > 2) resonances that can be excited by the I = 2 tide. 

With progressing nuclear evolution and consequent con- 
traction of the stellar core the frequency of the g-mode oscilla- 
tions increases, so that the 'gap' between potential resonances 
with strong retrograde g-modes and the strong r-modes widens. 
Thereby the retrograde strong g-modes disappear from the 
'tidal window' (a rcs > -mfl B ) so that spin-down of fast spin- 
ning early type binary stars depends on excitation of strongly 
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damped high radial order g-modes. This seems a not very effi- 
cient mechanism unless viscous effects (e.g. in the fi a -gradient 
zone outside the convective core) are significant or perhaps 
driving by inertial resonances in the convective core occurs. 
But, when in an eccentric binary system, the early type star is 
nearly corotating at periastron r-mode resonances are very ef- 
ficient in spinning the star further down to (pseudo)corotation. 
We expect a balance between spin-down by r-modes and spin- 
up by prograde g-modes excited by the high frequency compo- 
nents of the tidal forces near periastron. We intend to study this 
in a following paper, using the torque values determined here. 
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Appendix A: Tables 

Table Al. Fitting parameters for g^-mode resonances with an I = 2, m = tidal potential in a 10 M Q star with X c = 0.4. All 
tabulated values for the torque integral have been obtained for a fixed orbital separation. 

n s = o.o n c n s = 0.2 n c n s = 0.4 n c 
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Table A2. Fitting parameters for retrograde g^-mode resonances with an I = 2, m = 2 tidal potential in a 10 M Q star with 
X c = 0.4. All tabulated values for the torque integral have been obtained for a fixed orbital separation. 
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X 


10 42 


4 


.57 


X 


10" 


-4 


-0 


,8592 


-1 


,40 


X 


10 42 


5. 


.35 


X 


10" 


-4 


-0 


,9189 


-1 


.22 


X 


10 42 


4 


.96 


X 


10 


-4 


8 


-0 


.7501 


-1 


.93 


X 


1Q 40 


3 


.07 


X 


10" 


-3 


-0 


,8065 


-1 


.17 


X 


1Q 40 


3 


.51 


X 


10" 


-3 


-0 


,8693 


-8. 


.35 


X 


10 39 


3 


84 


X 


10 


-3 


9 


-0 


.7277 


-7. 


.76 


X 


10 42 


9 


.15 


X 


10" 


-5 


-0 


,7794 


-7, 


.71 


X 


10 42 


7 


48 


X 


10" 


-5 


-0. 


,8395 


-7. 


.92 


X 


10 41 


4, 


.30 


X 


10 


-4 


10 


-0 


.6434 


-2 


.01 


X 


10 43 


1. 


.92 


X 


10" 


-5 


-0 


,6985 


-9 


.26 


X 


10 42 


3 


.33 


X 


10" 


-5 


-0. 


,7602 


1. 


.27 


X 


10 42 


1 


.97 


X 


10 


-4 


11 


-0 


.5827 


-1 


,88 


X 


10 42 


5 


.97 


X 


10" 


-5 


-0. 


,6417 


-7, 


,80 


X 


10 41 


1 


.05 


X 


10" 


-4 


-0 


.7054 


-4 


.40 


X 


10 41 


1 


.25 


X 


10 


-4 


12 


-0 


.5525 


-3. 


.89 


X 


10 41 


2 


.51 


X 


10" 


-4 


-0. 


,6143 


-3 


,04 


X 


10 41 


2 


.75 


X 


10" 


-4 


-0 


.6795 


-1 


.30 


X 


10 41 


4 


.90 


X 


10 


-4 


13 


-0 


.5238 


-8, 


.97 


X 


1Q 40 


5 


.71 


X 


10" 


-4 


-0. 


,5849 


-2 


,58 


X 


10 41 


1. 


31 


X 


10" 


-4 


-0 


.6491 


-5 


.72 


X 


10 40 


4, 


.04 


X 


10 


-4 


14 


-0 


.4899 


-1. 


.29 


X 


10 41 


2 


.63 


X 


10" 


-4 


-0 


,5525 


-1 


.73 


X 


10 41 


1 


60 


X 


10" 


-4 


-0. 


,6170 


-9 


.25 


X 


10 40 


1 


.90 


X 


10 


-4 


15 


-0 


.4612 


-1, 


.47 


X 


1Q 40 


1. 


.44 


X 


10" 


-4 


-0 


5253 


-6 


.51 


X 


10 39 


2 


64 


X 


10" 


-4 


-0 


,5900 


-3 


.05 


X 


10 39 


2 


.33 


X 


10 


-4 


16 


-0 


.4390 


-7, 


.70 


X 


10 39 


4 


.13 


X 


10" 


-4 


-0 


,5046 


-6 


.79 


X 


10 39 


4. 


52 


X 


10" 


-4 


-0. 


,5695 


-5 


.67 


X 


1Q 39 


4, 


,15 


X 


10 


-4 


17 


-0 


.4259 


-2. 


.04 


X 


10 39 


6 


.86 


X 


10" 


-4 


-0 


,4919 


-1 


,24 


X 


10 39 


8. 


12 


X 


10" 


-4 


-0 


,5564 


-8 


.61 


X 


10 38 


9 


.36 


X 


10 


-4 


18 


-0 


.4101 


-1 


.02 


X 


10 38 


6 


.33 


X 


10" 


-4 


-0. 


,4757 


-1 


.15 


X 


10 38 


8. 


.91 


X 


10" 


-4 


-0 


5391 


-1 


.80 


X 


10 38 


6 


.93 


X 


10 


-4 


19 


-0 


.3935 


-6 


.68 


X 


10 38 


9 


.02 


X 


10" 


-4 


-0 


,4592 


-6 


.61 


X 


10 38 


6 


.68 


X 


10" 


-4 


-0 


,5216 


-3 


.90 


X 


10 38 


6 


.60 


X 


10 


-4 


20 


-0 


.3787 


-4, 


.97 


X 


10 37 


2 


.03 


X 


10" 


-3 


-0 


,4445 


-7, 


.44 


X 


10 37 


1 


.29 


X 


10" 


-3 


-0 


,5060 


-1 


.57 


X 


10 38 


9 


52 


X 


10 


-4 
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Table A3. Fitting parameters for prograde g|-mode resonances with an I = 2, m = 2 tidal potential in a 10 M Q star with 
X c = 0.4. All tabulated values for the torque integral have been obtained for a fixed orbital separation. 

n s = 0.2 n c n s = 0.3 n c n s = 0.4 n c 



k (T /Q c Tjrnfi (erg) Acr/Q. c go/Oc Tlm.o (erg) Aa/tt c up/flc Tim.fi (erg) A(t/O c 



1 


2 


1724 


1 


.92 


X 


10 47 


1 


.55 


X 


10" 


-6 


2 


.1519 


2 


,38 


X 


10 47 


1 


.43 


X 


10" 


-6 


2. 


.1330 


2. 


.98 


X 


10 47 


1 


26 


X 


10" 


-fi 


2 


1. 


.8006 


2 


.31 


X 


10 47 


1. 


.39 


X 


10" 


-6 


1. 


.7781 


1. 


.49 


X 


10 47 


3 


.53 


X 


10" 


-6 


1 


,7587 


6 


.71 


X 


10 46 


5 


.67 


X 


10" 


-6 


3 


1. 


.2773 


8. 


.31 


X 


10 46 


1. 


.11 


X 


10" 


-6 


1. 


.2623 


6 


.55 


X 


10 46 


1. 


.61 


X 


10" 


-6 


1 


,2498 


5 


.15 


X 


10 45 


2 


.28 


X 


10" 


-5 


4 





.9780 


7. 


.33 


X 


10 44 


2 


.94 


X 


10" 


-5 





.9638 


1. 


.37 


X 


10 45 


1. 


.87 


X 


10" 


-5 





,9532 


2 


.65 


X 


10 45 


1. 


.13 


X 


10" 


-5 


5 





.9011 


3 


.11 


X 


10 43 


1 


.33 


X 


10" 


-4 





.8836 


3 


.26 


X 


10 43 


1. 


.23 


X 


10" 


-4 





,8692 


6 


.30 


X 


10 43 


7 


.10 


X 


10" 


-5 


6 





.7629 


1. 


,46 


X 


10 44 


3 


.29 


X 


10" 


-5 





.7508 


2 


.57 


X 


10 44 


2 


.17 


X 


10" 


-5 


0. 


,7425 


2 


.56 


X 


10 44 


2 


.56 


X 


10" 


-5 


7 





.6783 


9 


.16 


X 


10 42 


1. 


.27 


X 


10" 


-4 





.6637 


1. 


.61 


X 


10 43 


8 


.03 


X 


10" 


-5 


0. 


,6534 


2 


.18 


X 


10 43 


6 


.48 


X 


10" 


-5 


8 





.6171 


1 


.00 


X 


10 42 


5 


.41 


X 


10" 


-4 





.6044 


3 


.29 


X 


10 42 


2 


.64 


X 


10" 


-4 





,5959 


6 


.81 


X 


10 42 


1. 


.62 


X 


10" 


-4 


9 





.6078 


2 


.55 


X 


10 42 


2 


.76 


X 


10" 


-4 





.5948 


1. 


.66 


X 


10 42 


2 


.75 


X 


10" 


-4 





,5853 


7 


.28 


X 


10 41 


3 


.55 


X 


10" 


-4 


10 





.5242 


5 


.54 


X 


10 43 


1. 


.11 


X 


10" 


-5 





.5152 


3 


.09 


X 


10 43 


2 


.21 


X 


10" 


-5 





,5103 


1 


.05 


X 


10 43 


7 


.07 


X 


10" 


-5 


11 





.4596 


9 


.80 


X 


10 42 


2 


.08 


X 


10" 


-5 





.4517 


3 


.54 


X 


10 42 


6 


.27 


X 


10" 


-5 





,4476 


4 


.44 


X 


10 42 


5 


.33 


X 


10" 


-5 


12 





.4213 


2 


.92 


X 


10 42 


4. 


.04 


X 


10" 


-5 





.4125 


5 


.07 


X 


10 42 


2 


.35 


X 


10" 


-5 





,4077 


6 


.97 


X 


10 42 


1. 


.70 


X 


10" 


-5 


13 





.3971 


6 


.25 


X 


10 41 


5 


.02 


X 


10" 


-5 





.3896 


2 


.15 


X 


10 41 


1. 


.19 


X 


10" 


-4 





,3855 


2 


.86 


X 


10 41 


6 


.81 


X 


10" 


-5 


14 





.3627 


4. 


.18 


X 


10 41 


6 


.77 


X 


10" 


-5 





.3565 


2 


.87 


X 


10 41 


8. 


.87 


X 


10" 


-5 





,3536 


1 


.76 


X 


10 41 


1. 


.30 


X 


10" 


-4 


15 





.3313 


3 


.67 


X 


10 39 


1. 


.56 


X 


10" 


-4 





.3256 


1. 


,99 


X 


10 39 


1. 


.69 


X 


10" 


-4 





,3232 


1 


.54 


X 


10 39 


1. 


.24 


X 


10" 


-4 


16 





.3057 


5 


.00 


X 


10 39 


2 


.52 


X 


10" 


-4 





.3004 


3 


.85 


X 


10 39 


2 


.43 


X 


10" 


-4 





,2981 


3 


.26 


X 


10 39 


2 


.33 


X 


10" 


-4 


17 





.2896 


9 


.20 


X 


10 39 


2 


.98 


X 


10" 


-4 





.2840 


9 


.48 


X 


10 39 


2 


.83 


X 


10" 


-4 





,2813 


9 


.43 


X 


10 39 


2 


.65 


X 


10" 


-4 


18 





.2758 


2 


,06 


X 


10 38 


6 


.43 


X 


10" 


-4 





.2711 


3 


.34 


X 


10 38 


6 


.40 


X 


10" 


-4 


0. 


,2690 


6 


.06 


X 


10 38 


5 


.35 


X 


10" 


-4 


19 





.2586 


5 


.45 


X 


10 38 


1. 


.08 


X 


10" 


-3 





.2545 


5 


.09 


X 


10 38 


1. 


.01 


X 


10" 


-3 


0. 


,2529 


4, 


,88 


X 


10 38 


9 


.27 


X 


10" 


-4 


20 





.2421 


6 


,70 


X 


10 37 


1 


.66 


X 


10" 


-3 





.2384 


8. 


.80 


X 


10 37 


1 


.47 


X 


10" 


-3 





,2370 


1 


.05 


X 


10 38 


1 


.28 


X 


10" 


-3 



Table A4. Fitting parameters for r|-mode resonances with an I = 2, m = 2 tidal potential in a 10 M Q star with X c = 0.4. All 
tabulated values for the torque integral have been obtained for a fixed orbital separation. 



n s = 0.1O c O s = 0.2f2 c 
k go/Oc Tim,o (erg) Aa/Q c go/Op 7I m , (erg) Aa/Q c 






-3 


.3237 


X 


10" 


-2 


-5.63 


X 


10 41 


5.71 


X 


10" 


■6 


-6 


.6098 


X 


10" 


-2 


-2 


.57 


X 


10 41 


7.35 


X 


10" 


-5 


1 


-3 


.3191 


X 


10" 


-2 


-7.69 


X 


10 40 


4.80 


X 


10" 


■5 


-6 


.5845 


X 


10" 


-2 


-1 


.20 


X 


10 41 


2.40 


X 


10" 


-4 


2 


-3 


.3133 


X 


10" 


-2 


-1.04 


X 


10 41 


3.61 


X 


10" 


-5 


-6 


,5292 


X 


10" 


-2 


-1 


.41 


X 


10 41 


2.03 


X 


10" 


-4 


3 


-3 


.3031 


X 


10" 


-2 


-8.46 


X 


10 40 


3.03 


X 


10" 


-5 


-6 


.4479 


X 


10" 


-2 


-1 


.03 


X 


10 41 


1.87 


X 


10" 


-4 


4 


-3 


.2896 


X 


10" 


-2 


-4.67 


X 


10 40 


3.14 


X 


10" 


■5 


-6 


.3467 


X 


10" 


-2 


-5 


.42 


X 


10 40 


1.85 


X 


10" 


-4 


5 


-3 


.2730 


X 


10" 


-2 


-2.43 


X 


10 40 


3.46 


X 


10" 


-5 


-6 


.2277 


X 


10" 


2 


-2 


.96 


X 


10 40 


2.00 


X 


10" 


-4 


6 


-3 


.2534 


X 


10" 


-2 


-1.27 


X 


10 40 


3.97 


X 


10" 


-5 


-6 


.0922 


X 


10" 


-2 


-1 


.62 


X 


10 40 


2.11 


X 


10" 


-4 


7 


-3 


.2315 


X 


10" 


-2 


-7.93 


X 


10 39 


4.26 


X 


10" 


-5 


-5 


.9485 


X 


10" 


-2 


-9 


.14 


X 


10 39 


2.38 


X 


10" 


-4 


8 


-3 


.2084 


X 


10" 


-2 


-4.02 


X 


10 39 


4.63 


X 


10" 


■5 


-5 


.8029 


X 


10" 


-2 


-5 


.19 


X 


10 39 


2.51 


X 


10" 


-4 


9 


-3 


.1842 


X 


10" 


-2 


-1.48 


X 


10 39 


5.51 


X 


10" 


-5 


-5 


,6565 


X 


10" 


-2 


-1 


,89 


X 


10 39 


3.12 


X 


10" 


-4 


10 


-3 


.1581 


X 


10" 


-2 


-2.24 


X 


10 38 


8.18 


X 


10" 


-5 


-5 


.5121 


X 


10" 


-2 


-3 


.39 


X 


10 38 


7.66 


X 


10" 


-4 



fi s =0.3 fi c Q s = 0.4fi c 
k op/^lz Trnfi (erg) Aa/Q c ao/Q c ?Im,o (erg) Aff/ii c 






-9.8117 


X 


10" 


-2 


-1.74 


X 


10 41 


3.87 


X 


10" 


-4 


-1.2898 


X 


10" 


-i 


-3.26 


X 


10 41 


4.25 


X 


10" 


-4 


1 


-9.7284 


X 


10" 


-2 


-1.12 


X 


10 41 


9.12 


X 


10" 


-4 


-1.2685 


X 


10" 


-l 


-2.23 


X 


10 41 


6.54 


X 


10" 


-4 


2 


-9.5530 


X 


10" 


-2 


-1.95 


X 


10 41 


3.78 


X 


10" 


-4 


-1.2328 


X 


10" 


-l 


-4.60 


X 


10 41 


3.18 


X 


10" 


-4 


3 


-9.3144 


X 


10" 


-2 


-1.51 


X 


10 41 


3.93 


X 


10" 


-4 


-1.1833 


X 


10" 


-l 


-2.73 


X 


10 41 


4.26 


X 


10" 


-4 


4 


-9.0057 


X 


10" 


-2 


-4.14 


X 


10 40 


7.76 


X 


10" 


-4 


-1.1145 


X 


10" 


-l 


-3.08 


X 


10 40 


2.42 


X 


10" 


-3 


5 


-8.6610 


X 


10" 


-2 


-2.93 


X 


10 40 


4.95 


X 


10" 


-4 


-1.0515 


X 


10" 


-l 


-4.60 


X 


10 40 


5.49 


X 


10" 


-4 


6 


-8.2982 


X 


10" 


-2 


-2.23 


X 


1Q 40 


4.14 


X 


10" 


-4 


-9.8842 


X 


10" 


-2 


-3.21 


X 


10 40 


5.26 


X 


10" 


-4 


7 


-7.9327 


X 


10" 


-2 


-9.10 


X 


10 39 


5.65 


X 


10" 


-4 


-9.2639 


X 


10" 


2 


-7.71 


X 


10 39 


1.12 


X 


10" 


-3 


8 


-7.5855 


X 


10" 


-2 


-5.14 


X 


10 39 


6.24 


X 


10" 


-4 


-8.6945 


X 


10" 


2 


-5.05 


X 


10 39 


1.08 


X 


10" 


-3 


9 


-7.2437 


X 


10" 


-2 


-2.03 


X 


10 39 


7.43 


X 


10" 


-4 


-8.1781 


X 


10" 


-2 


-3.74 


X 


10 39 


7.49 


X 


10" 


-4 


10 


-6.9712 


X 


10" 


-2 


-2.28 


X 


10 38 


3.98 


X 


10" 


-3 


-7.7343 


X 


10" 


-2 


-3.83 


X 


10 38 


1.02 


X 


10" 


-3 



